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ABSTRACT 

Compton scattering of X-rays in the bulk flow of the accretion column in magnetic cat- 
aclysmic variables (mCVs) can significantly shift photon energies. We present Monte 
Carlo simulations based on a nonlinear algorithm demonstrating the effects of Comp- 
ton scattering on the H-like, He-like and neutral Fe Ka lines produced in the post- 
shock region of the accretion column. The peak line emissivities of the photons in 
the post-shock flow are taken into consideration and frequency shifts due to Doppler 
effects are also included. We find that line profiles are most distorted by Compton 
scattering effects in strongly magnetized mCVs with a low white dwarf mass and high 
mass accretion rate and which are viewed at an oblique angle with respect to the ac- 
cretion column. The resulting line profiles are most sensitive to the inclination angle. 
We have also explored the effects of modifying the accretion column width and using 
a realistic emissivity profile. We find that these do not have a significant overall effect 
on the resulting line profiles. A comparison of our simulated line spectra with high 
resolution C/iandra/HETGS observations of the mCV GK Per indicates that a wing 
feature redward of the 6.4 keV line may result from Compton recoil near the base of 
the accretion column. 

Key words: accretion - line: profiles - scattering - binaries: close - white dwarfs - 
X-rays 



1 INTRODUCTION where /j, is the mean molecular weight and a; s is the 

. , , , . , . shock height. For WD masses M wd » (0.5 - 1.0) Mp, 

Magnetic cataclysmic variables mCVs) are close interact- , , . , ™ r , , , , . , r 

... r /,,,t-n\ anc i typical mCV parameters, the shock temperature is 

mg binaries consisting ol a magnetic white dwarf (WD) and IrTn /,„ ,„,, /, , s ir mi i 

f a a rFTr Jinnch w v.- a kTs ~ ( 10 _ 40) keV ps (1 - 4) x 10 K. The plasma 



a low mass red dwarf IjWarnerll 19951 ). Near the white dwarf 
surface the accretion flow in mCVs is confined by the mag- 



in the post-shock region of the column cools by emitting 
bremsstrahlung X-ray s and optical/IR cy clotron radiation 
netic field of the WD and is channelled to the magnetic j Lamb % Mastcrs gj| | K in g fc Lasotalfl979ll . Since cooling 
pole region s) ot the WD. torming an accretion column sec , „ , , , . . , .„ , . 

i ; ~ — t l „„ J II ' 'j i „ T n, i rrz n i„ J „ occurs along the how, the post-shock region is stratified m 

gOB l2™dS; @H; \Wuet^&M, for density a J\ emperature . £ he heigM of S the pogt . shock re _ 

gion is determined by the cooling length. For a flow with 
Near the base of the accretion column, material m su- , , , , ,. ,, , 

. , „' _ . only brem sstrahlung cooling, the shock height is given by 

persomc free-fall is brought to rest on the WD surface, form- j W u. Chanmugam fe ShavivS 
mg a standing shock which heats and ionizes the accreting 

plasma. The shock temperature T s depends mainly on the / ^ \ f M d \ ^ 

mass A/ wd and radius i? wd of the WD and is given by (e.g. x * ~ 3 x 10 _ 2 _i TmH 

rw3H) Vlgcm s A ®/ 2 

fcTs= 3GAWm H ( (1) x (l^) cm ( 2 ) 

where rh is the specific mass accretion rate. 

A plasma temperature of kT ~ 10 keV is sufficient to 
* E-mail:aimee, kuncic@physics.usyd.edu.au fully ionize elements such as argon, silicon, sulphur, alu- 

f Centenary Fellow minium or calcium. Heavier elements such as iron can be 
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highly ionized, resulting in H-like Fe xxvi and He-like Fe 
xxv ions. K-shell transitions in Fe xxvi and Fe xxv ions 
give rise to Ka lines at 6.97 keV and 6.6 75 keV, respectively 
(see e.g. I Wu. Cropper fe Ram say 2001). The irradiation of 
low ionized and neutral iron by X-rays above the Fe K edge 
produces fluorescent Ka emission at approximately 6.4 keV 
inside the WD atmosphere, beneath the accretion column 
and in surrounding areas. The natural widths of the Fe Ka 
lines are small, but the lines can be Doppler broadened by 
the bulk and thermal motions of the emitters in the post- 
shock flow. The bulk velocity immediately downstream of 
the shock is ^ 0.25 (GM wd /i? wd ) 1/2 ~ 1000 kms" 1 for typ- 
ical mCV parameters. Lines can also be broadened by scat- 
tering processes. Compton (electron) scattering is expected 
to be more impo rtant than resonance (ion) scattering f or the 
Ka transitions l|Pozdnvakov. Sobol fe Sunvaevl fl983). For 
an mCV with specific accretion rate m ~ 10 gcm -2 s _1 
the electron number density is n c ~ 10 16 cm -3 for a shock 
heated region of thickness x B ~ 10 7 cm, giving a Thompson 
optical depth of t ~ 0.1. Thus, one in every ten photons 
would encounter an electron before leaving the post-shock 
region. The relative importance of Doppler shifts, thermal 
Doppler broadening and Compton scattering depends on the 
ionization structure in the post-shock flow. 

X-ray observations by Chandra /HETGS and 
ASCA/S1S have reveale d significant broadening of some 
Fe Ka lines in mCV s (|Hellier. Mukai. fe Osborne.! Il998t 
iHellier fe Mukail | 20041 ) . It was suggested that Compton 
scattering in the accretion column is largely responsible for 
the broadening in the observed lines. Doppler broadening 
should only be significant in lines emitted close to the 
shock. The absence of Doppler shifts in the observed H-like 
and He-like lines suggests that these photons may be 
emitted predominantly from region s of lower velocity nea r 
the base of the accretion column (|Hellier fe Mu kai 2004). 
The observed line profiles have yet to be fully interpreted 
with a quantitative model that takes into account Compton 
scattering effects in a complex ionization structure. 

In this paper, we study the effects of Compton scat- 
tering in the accretion colum n of mCVs u sing a non- 
linear Monte Carlo algorithm l|Cullenl |20bTal lbh that self- 
consistently takes into account the density, velocity and tem- 
perature structure in the column (|Wu et al.ll200ll ). The ef- 
fects of dynamical Compton scattering are also included . In 
a preliminary investigation (|Kuncic. Wu fe C ullcn 2005|), it 
was found that Fe Ka line photons emitted from the dense 
base of the accretion column undergo multiple Compton 
scatterings and as a result, the base of the line profile is 
substantially broadened. Photons emitted near the shock 
can also undergo scatterings with hot electrons immediately 
downstream of the shock, as well as cold electrons in the 
pre-shock flow before escaping the column. The resulting 
line profiles display a shoulder-like feature redward of the 
line centre. More significant broadening is observed when cy- 
clotron cooling is sufficientl y strong t o produce a dense, com- 
pact post-shock region (see lWull2000l . for example). Here, we 
make three substantial improvements to the previous study: 
(i) we include Doppler effects; (ii) the photon source regions 
in the post-shock column are determined from the ioniza- 
tion structure, rather than specified arbitrarily; and (iii) the 
effects of different viewing angles are fully explored. The 
paper is organized as follows: the theoretical outline and ge- 



ometry of the model are described in Section 2. Numerical 
results for cases where cyclotron cooling is negligible and 
when it dominates are presented and discussed in Section 3. 
A summary and conclusions are presented in Section 4. 



2 THEORETICAL MODEL 
2.1 Physical Processes 

Line photons can undergo energy changes when scattering 
with electrons, resulting in distortions in the line profile. The 
ene rgy change of a photon per scattering is given by (e.g. 
see iPozdnvakov et al.|[l983h . 
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where E is the initial photon energy, 7m e c 2 is the electron 
energy, with 7 = (1 — /3 2 ) -1 / 2 and where /3 = v/c includes 
both thermal and bulk motion, 8 — cos" 1 fi is the incident 
photon propagation angle measured relative to the electron's 
direction of motion, and the scattering angle is a. The prime 
superscript denotes quantities after a scattering event. Al- 
though the energy change per scattering is typically small, 
the line profile can be broadened considerably as a result of 
multiple scatterings if the optical depth is large. 

Photons scattering with hot electrons (i.e. kT c > E c 
where T e is the electron temperature and E c is the line centre 
energy) will gain energy, while photons scattering with cold 
electrons (kT e <C E c ) will lose energy due to recoils. In the 
mCV context, Compton recoil can be important for photons 
scattering with electrons in the cold pre-shock flow and also 
near the base of the column, where f3 rapidly decreases and 
where the optical depth is high. From equation © , the frac- 
tional energy change due to Compton recoil (E ^> 7m e c 2 ) 
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(1 — cos a) 



(4) 



In the post-shock accretion column in mCVs, line pho- 
tons undergo thermal Doppler broadening as well as Doppler 
shifts. The bulk velocity of the accreting material immedi- 
ately downstream of the shock is ~ 1000 kms -1 . For all 
inclination angles i (see Figure [IJ, the bulk flow is moving 
away from our line of sight, so the line centre energy is red- 
shifted by an amount AE/E ~ /3cosi. However, since the 
bulk velocity in the post-shock flow in mCVs is always less 
than a few x 1000 kms -1 , giving j3 < 10 -2 , Doppler shifts 
are expected to be negligible. Thermal Doppler broaden- 
ing, on the other hand, is expected to be of order AE/E ~ 
yj (2kT B /rrnc 2 ) ~ 0.002 for lines emitted in the hottest re- 
gions of the post-shock flow (i.e. immediately downstream 
of the shock), where mi is the mass of the ion. 

In strongly magnetized mCVs, cyclotron emission is the 
dominant cooling process. The effect of this additional elec- 
tron cooling process in the post-shock region is to reduce 
the shock height and m odify the densi ty and temperature 
structure of the region (|Wu et alJll994l ). This can enhance 
Com pton scattering features in line profiles (|Kuncic et al.l 
l2005h . 



Compton Scattering of Fe Ka Lines in mCVs 3 



Accretion Flow 




Figure 1. Schematic illustration of the geometry of the magne- 
tized white dwarf accretion column showing approximate loca- 
tions of Fe Ka source regions. 

2.2 Geometry of the Accretion Column 

The geometry of an mCV accretion column is shown in fig- 
ure [T] The column is modelled as a cylinder and is divided 
into a shock heated region and a cool pre-shock region. In 
the pre-shock region, the density and velocity of the accret- 
ing material are constant. The velocity of the pre-shock flow 
is approximated by the free-fall velocity at the shock height, 
vs(x s ) = [GM w d/(l?wd + i s )]^ 2 , and the plasma is assumed 
to be cold (kT a <C 1 keV). The electron number density in 
the pre-shock flow is n c = m/((iraHi)i), where rh is the spe- 
cific mass accretion rate. In the post-shock region, the den- 
sity, velocity and temperature profiles ar e calculated usin g 
the hydrodynamic solution described in IWu et all (| 1994T ) ■ 
for bremsstrahlung and cyclotron cooling. The parameters 
that determine the structure of the post-shock region are 
the WD mass M w d, the WD radius i? w d, the specific mass 
accretion rate rh and the ratio of the efficiencies of cyclotron 
cooling to bremsstrahlung cooling e s . 

The accretion column is viewed from an inclination an- 
gle i, measured relative to the column axis (see Fig. [TJ. An 
inclination angle i — 0° corresponds to viewing the column 
along its axis, towards the WD, while viewing the column 
from an inclination angle i — 90° is equivalent to viewing 
the column from the side. 

The line photons are injected into the post-shock region 
at a specific dimensionless height ( = x/x s at or above the 
WD surface, where f = and £ = 1 correspond to the 
WD surface and shock surface, respectively (see Fig.[T]). For 
each set of mCV parameters, the injection point of the line 
corresponds to the location in the post-shock region where 
the emissivity of the specific line p eaks, acco r ding to the 
ionization structure determined by IWu et all lj200lf ). The 
bulk and thermal velocities at the line injection height £ are 



used to calculate the Doppler shift and broadening of the 
lines. 

The Monte Carlo technique is used to model Comp- 
ton scattering effects on photon propagation in the column. 
The distance to a tentative scattering point is determined 
using an algorithm ba sed on a nonlinear transport tech- 
nique (|Stern et al.lll995h which integrates t he mean free p ath 
over the spatially varying electron density (|Cullenll2001a| [bT). 
The scattering cross-section is determined from the Klein- 
Nishina formula and the momentum vector at the scattering 
point is drawn from an isotropic Maxwellian distribution at 
the local temperature. A rejection algor ithm is used t o de- 
cide whether the scattering is accepted l|Cullenll2001bh . For 
an accepted event the energy and mom entum changes of 
the photon are calculated as described in lPozdnvakov et al.l 
(1983). In each simulation photons are followed until they 
leave the column and binned to form a spectrum. A full de- 
scriptio n of the numerical algorithm can be found in ICullenl 
1(2001 tf) . 

In the simulations described below, we make the follow- 
ing simplifying assumptions: a fixed number of 10 s photons 
are used to simulate each line; a single energy is used for 
each line (whereas in reality, the neutral and Lymann a 
transitions are doublets, with energies 6.391/6.404 keV and 
6.952/6.973 keV, and the He- like transition has both reso- 
nant, inter-combination and forbidden components); a fixed 
accretion column width is used and a single line injection 
site is used. 



3 RESULTS AND DISCUSSION 

We present simulations of Compton scattering of Fe Ka 
lines in an mCV accretion column for two different WD 
mass-radius values: M w d = 0.5 Mq, -R w d = 9.2 x 10 s cm 
and M wd = 1.0 Mq, R vd = 5.5 x 10 s cm (|Nauenberg||l972l ). 
For each mass, we consider two different specific mass ac- 
cretion rates, m = 1 gcm _2 s _1 and m = 10 gcm _2 s _1 , 
and inclination angles i = 0°, 45° and 90° with ±5° range. 
The cross-sectional radius of the accretion column is fixed 
at 0.1-R w d. We investigate the effect of varying the column 
width in section 3.3. 

Figure [2] shows the temperature, velocity and density 
profiles of an accretion column for an mCV with WD mass 
M w d = 0.5 Mq and rh = 10 gem -2 s _1 , for cases where the 
ratio of cyclotron to bremsstrahlung cooling at the shock is 
e s = 0, 10, and 100. The temperature decreases monotoni- 
cally from a shock temperature kT s ~ 14keV (£ = 1), to a 
small finite value at the base of the column (£ = 0).The 
mass density in the post-shock region is determined by 
P(C) = 4m/«b(C) an d is a minimum at the shock and reaches 
a maximum at the base of the column. The bulk velocity of 
the accreting material at the shock is 0.25i;ff and decreases 
to zero at the base of the column, where the plasma settles 
on the WD surface. 

Figures [3] and [4] show the simulated line spectra 
for the case where cyclotron cooling is negligible and 
bremsstrahlung cooling dominates (e s = 0). Figure [5] shows 
the simulated profiles for M w d = 0.5 Mq and M w d = 
1.0 Mq when cyclotron cooling dominates bremsstrahlung 
(e s = 10). In these cases, the neutral Fe Ka line is emitted 
at the WD surface (£ = 0), where the bulk velocity is ap- 
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Figure 2. The profiles of (a) electron temperature T e , (b) mass density p e and (c) bulk velocity v-^ = p^c in the post-shock accretion 
column of an mCV with M wd = 0.5Mq and fa = 10gcm~ 2 s^ 1 . £ = corresponds to the base of the accretion column and £ = 1 to 
the shock. The solid line shows the profile for the case where the ratio of cyclotron to bremsstrahlung cooling at the shock is e B = 0, the 
dotted line for e s = 10 and the dashed line for e B = 100. In (a) the corresponding mean thermal electron velocity /3 e = v c /c is also shown 
as a function of 



proximately zero and the thermal velocity of the plasma is 
small. The 6.675 keV Fe Ka line is emitted from the lowest 
few percent of the column (£ ~ 0.003) where the velocity 
of the infalling material and the thermal electron velocity 
are still relatively small. These lines thus show very little 
Doppler broadening. The 6.97 keV line, however, is emit- 
ted much closer to the shock (£ ~ 0.16), in regions where 
the temperature of the accreting material is considerably 
hig her and thus suff ers more substantial Doppler broadening 
(see lWu et al.ll200ll ). The simulated line spectra for the case 
where the accretion column radius is fixed at 4.6 x 10 7 cm 
is shown in Figure [S] Figure [7] shows the simulated pro- 
files for the case where the line photon injection along the 
flow is specified according to the emissivity profile model of 
IWu et all (|200ll ). In figure [5] we compare our simulated line 
spectra with an observation of the mCV GK Per detected 
by Chandra /HETGS. 

Table [T]shows a comparison of the FWHM of the Fc Ka 
lines for M wd = 1.0 Mq and M wd = 0.5 Mq with low and 
high fn and for e B = and e s — 10. 



3.1 No Cyclotron Cooling 

When cyclotron cooling is negligible (e s = 0), the 6.4, 6.675 
and 6.97 keV lines are mostly emitted at heights above the 
WD surface of x » 0, and 1.3 x 10 6 cm for M wd = 1.0 Mq 
and i«0, 1 .0 X 10 5 , and 5. 1 x 10 6 cm for M wd = 0.5 Mq, 
respectively (|Wu et al.ll200ll) . Note that Doppler broaden- 
ing is most significant for the 6.97 keV line which is emit- 
ted further up in warmer regions of the post-shock column. 
Figures [3] and [3] show that in all cases, there is substantial 
Compton broadening near the base of the line profiles due to 
scatterings in the post-shock region where photons may lose 
and gain energy. For example, for the 6.4 keV line emerging 
at i = 90° in Fig. [3^, the Compton broadened wings con- 
tain « 28% of the total line photons. The blue wing extends 
up to AE » 0.09 keV above the line centroid, while the 
red wing extends down to AE » 0.16 keV (see Appendix). 
Compton features are generally more prominent at higher 
accretion rates since this gives higher optical depths in the 
column and hence, a higher scattering probability. For com- 
parison, the average number of scatterings per photon in the 



Table 1. FWHM (eV) values of the (a) 6.4 keV, (b) 6.675 keV 
and (c) 6.97 keV lines calculated for two different values of fa and 
for M wd = 1.0 Mq and M wd = 0.5 Mq. The ratios of cyclotron 
to bremsstrahlung cooling at the shock are e B = and e s = 10. 



(a) 
6.4 keV 


M wd = l.OM 
fa = 1 (10) gem -2 s _1 


M wd = O.5M 
fn = 1 (10) gem -2 s _1 


e B = 
e a = 10 


5(7) 
4(5) 


4(5) 
3(4) 




(b) 
6.675 keV 


M wd = 1.0 Mq 

fa = 1 (10) gem -2 s~ x 


M wd =0.5 Mq 
fa = 1 (10) gem - s — 1 


e B = 
e H = 10 


6(7) 
5(6) 


4(5) 
3(5) 




(c) 
6.97 keV 


M wd = 1.0 Mq 

fa = 1 (10) g cm -2 s _1 


M wd = 0.5 M 
fa = 1 (10) gem -2 s — 1 


e B = 
e H = 10 


8(9) 
7(15) 


10(10) 
8(10) 



fa = lgcm~ 2 s _1 case is 0.5 but this increases to 2 in the 
fn = 10 g cm -2 s _1 case. 

For M wd = 0.5 Mq (Fig. [4j , a recoil tail redward of the 
line centroid can be seen in the cases i = 45° and 90° specif- 
ically for high fn (Fig. |4p) The shock temperature for these 
cases is kT B ~ 14 keV, while for the M wd = 1.0 Mq case, the 
shock temperature is much higher (kT B » 34keV for low fn 
and kT B ~ 46 keV for high fn). The plasma temperature in 
the 1.0 Mq case is too high to produce any downscattering 
features and a large M wd (small 7? wd ) has an accretion col- 
umn with a lower optical depth, hence the absence of these 
features in Fig. [3] Again the downscattering features are 
more enhanced in the high fn cases since the optical depth 
is higher (c.f. Figs. [4ji and|4p for i = 90°). Note also that 
the 6.97 keV line centre is Doppler shifted slightly redward 
which is most noticeable for the 0.5 Mq case. 

In addition to the recoil signatures, upscattering fea- 
tures are also seen in the line spectra for the high fn cases 
(Figs. [3Jd and |4p) when i = 90°. For these inclination an- 
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Figure 3. Profiles of Fe Ka lines scattered by electrons in the accretion column of an mCV with M wd = 1.0 Mq and with (a) 
rh = 1 gcm -2 s — 1 and (b) rh = 10 gem - 2 s — 1 . The shock temperature kT B is indicated and the ratio of cyclotron to bremsstrahlung 
cooling is e s = 0. 



gles, photons emerge from the column with a final scattering 
angle cosine // ~ 0. Equation [3] then gives 

J « 1 - ■ (5) 

Since head on collisions (/x < 1) have a higher probability, 
E' ~ (l + \n\(3)E. The upscattering features are more promi- 
nent in Fig. [3Jd than in Fig. |4p because the shock temper- 
ature, and hence /3, is higher. The sharpness of the upscat- 
tering features seen in the 6.4 keV profiles in particular can 
be attributed to the maximum possible energy gain when 
fl w — 1 and to very few hot electrons with energies beyond 
a few standard deviations of the mean thermal (Maxwellian) 
energy. 

Scattering features are less prominent in the low m 
cases (Figs. [5^ and > since the overall optical depth is 
smaller. For Af w d = 1.0 Mq and rh — 1 gcm -2 s -1 for in- 
stance, the optical depth across the column at the shock 
is t ~ 0.04, which is a factor 10 smaller than that for the 
rh = 10 gcm~ 2 s -1 case (since r oc n c cx rh). For higher 
WD masses (i.e. smaller J? w d), the optical depth across the 
column is smaller, so fewer photons are scattered, especially 
for low rh cases. This explains the difference in scattering 
features in the profiles shown in Figs. [3] and [4] (in particular, 
the high rh, i = 90° cases). 

Photons observed at an inclination angle of i — 0° prop- 
agate through the entire length of the column before es- 
caping. These photons propagate through a thick section of 
cold pre-shock flow and can thus downscatter, resulting in 
broadening red ward of the line centre. However, the electron 



number density, and hence, optical depth, in the pre-shock 
flow is small so recoil effects are correspondingly small. Fur- 
thermore, fewer photons are detected at i — 0° because the 
solid angle centered around i — 0° is considerably smaller 
than that centered around larger inclination angles. 



3.2 Cyclotron Cooling Dominated Flows 

Figure [S] shows the profiles of Fe Ka lines emitted from 
mCVs with M wd = 1.0 Mq and M wd = 0.5 Mq where 
i — 90°, rh — 10 gcm _2 s _1 and the ratio of cyclotron to 
bremsstrahlung cooling at the shock is e s = 10. Also plotted 
is the thermal Doppler broadened line profile before scat- 
tering (dotted curve). The additional cooling of the plasma 
in the post-shock flow enhances the density and results in 
larger scattering optical depths. Because the electron tem- 
perature of the shock-heated region decreases as e s increases, 
the peak emissivities of the 6.675 keV and 6.97 keV Fe lines 
are found at larger values of £, as can be seen in Fig. [2] which 
shows the temperature, density and velocity profiles for dif- 
ferent cases of e a . The additional cyclotron cooling results in 
a decrease in the electron temperature and bulk velocity in 
the post-shock column at a fixed £. 

For e s = 10, the 6.4, 6.675 and 6.97 keV lines are emit- 
ted from heights above the WD surface of x » 0, 3.1 x 10 5 
and 2.4 x 10 7 cm for M wd = 0.5 Mq, and x w 0, and 

5.3 x 10 7 cm for M wd = 1.0 Mq, respectively (|Wu et all 
2001). The emission height of the 6.97 keV line is signifi- 
cantly higher than in the e s = case and consequently, 
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(a) iW wd = 0.5 Mq, rh = 1 gem 2 s 1 , kT s = 14keV, no cyclotron cooling 
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Figure 4. Profiles of Fe Ko lines scattered by electrons in the accretion column of an mCV with M wd = 0.5 Mq and with (a) 
rh = 1 gcm -2 s _1 (top panels) and (b) rh = 10 gcm _2 s" 1 (bottom panels). The shock temperature kT B is indicated and the ratio of 
cyclotron to bremsstrahlung cooling is e s = 0. 




Figure 5. Simulated profiles of Compton scattered Fe Ko lines in the accretion column of an mCV with (a) Af wd = 1.0 Mq and (b) 
M w j = 0.5 Mq. In both cases the inclination angle is i = 90°, the specific mass accretion rate is rh = 10 gem -2 s — 1 and the ratio of 
cyclotron to bremsstrahlung cooling at the shock is e s = 10. The dotted curves show the contribution to each line profile by thermal 
Doppler broadening only. 



thermal Doppler broadening is more prominent. The addi- 
tional cooling in the accretion column also results in a lower 
plasma temperature at the base of the column, £ = (see 
Fig. [2} , near where irradiation of neutral iron occurs and the 
fluorescent 6.4 keV line is emitted. This line thus undergoes 
less thermal Doppler broadening than in the e s = case. As 
a result of Compton scattering, the line profiles show addi- 
tional broadening and recoil tails redward of the line centre 
(Fig. [SJd) . Upscattering features are also evident in the 6.4 
and 6.675 keV lines especially for Af w( j = 1.0 Mq (Fig. [5^,). 

Overall, Compton scattering features are more promi- 



nent for cyclotron cooling dominated accretion columns, 
especially for high accretion rates. For photons viewed at 
i — 90°, upscattering and recoil features that are only 
marginally seen in the e s = case are considerably more 
pronounced in the e = 10 case. Thus, we expect Compton 
scattering features to be most conspicuous in Fe Ka lines 
emitted in strongly magnetized mCVs accreting at a high 
rate. 
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Figure 6. Simulated profiles of Compton scattered Fe Ka lines in an accretion column of an mCV with (a) M wd = 1.0 Mq and (b) 
M wd = 0.5 Mq both with a specific mass accretion rate of m = 10 gem -2 s —1 , a ratio of cyclotron to bremsstrahlung cooling of e s = 10, 
an inclination angle of i = 90° and accretion column radius of 4.6 X 10 7 cm . The dotted curves show the profiles when the column radius 
is fixed at 0.1-R wd (equivalent to the profiles shown in Fig[5jl. 



3.3 Effect of Accretion Column Radius 

The width of the accretion column in mCVs is poorly known 
and may vary significantly from system to system. In the 
results presented so far, we have fixed the accretion col- 
umn radius to be 0.1-R w d- Here, we investigate the effect 
of relaxing this assumption. Figure [6j shows the resulting 
profiles of Fe Ka lines emitted from two mCVs of different 
masses (M wd = 1.0 Mq and M wd = 0.5 Mq) but with the 
same absolute accretion column radius of 4.6 x 10 7 cm. This 
corresponds to 0.08i? w d for the 1.0 Mq case (Fig. [6^,) and 
0.05ii w d for the 0.5 Mq case (Fig. [Bp). The other parame- 
ters used to produce the profiles are the same as those used 
for Fig. \5\ namely m = 10 gem -2 s -1 , e s = 10 and i = 90°. 
Fig.[5Jshows that the changes in the accretion column width 
generally result in changes near the base of the line profiles, 
which is broadened by multiple scatterings. The effect of in- 
creasing (decreasing) the width of the accretion column is 
to increase (decrease) the overall scattering optical depth. 

The Thompson optical depth across the base of the ac- 
cretion column is ~ 20 for the 1.0 Mq case (Fig. [6ji) and 
~ 40 for the 0.5 Mq case (Fig. Ep)- In comparison, the 
Thompson optical depth across the column base for the 
same lines when the column has a width of 0.1i? w d (dot- 
ted curves in Fig. [5J are ~ 25 for M wd — 1.0 Mq and ~ 80 
for M w d = 0.5 Mq. Fig. [6j indicates that due to the nonlin- 
ear nature of multiple scatterings in the accretion column, 
particularly near the base, small changes in optical depth 
associated with the accretion column geometry can result 
in significant changes in line profiles. These changes mostly 
effect the base of the line profiles. 



3.4 Emissivity Profile Effects 

Realistically, the photon source regions are determined by 
an emissivity profile along the post-shock region of the ac- 
cretion column. In the results presented so far, all the pho- 
tons were injected at a single height in the accretion column 
corresponding to the locati on of the peak i n the emissiv- 
ity profile, as calculated by IWu et all 120011 ). Here, we in- 
vestigate the effect of spreading the photon injection site 
over a finite range of heights in the post-shock region ac- 
cording to the calculated emissivity profiles. Only the 6.675 



W 2 fT 
10* r 




in -6 1 

6.4 6.5 6.6 6.7 6.8 6.9 7.0 7.1 7.2 
Energy (keV) 

Figure 7. Simulated profiles of 6.675 keV and 6.97 keV Fe lines 
in an mCV accretion column with M wd = 1.0 Mq, a specific 
mass accretion rate of m = 1 gcm^ 2 s~ 1 , no cyclotron cooling 
and viewed at an inclination angle of i = 45° (solid lines). An 
emissivity profile is used to disperse the photon injection site 
over a finite range of heights along the post-shock column. The 
dotted curves show the corresponding profiles for photons injected 
at a single height where the emissivity peaks for the same mCV 
parameters (equivalent to the profiles shown in Fig |3ji, middle 
panel). 



keV and the 6.97 keV line are studied in this manner since 
much of the fluorescence 6.4 keV yield derives from beneath 
the column, where X-ray irradiation is strongest so the in- 
jection site is always at the base of the column. Figure [7] 
shows an example of the emissivity profile effect for the 
6.675 keV and 6.97 keV Fe Ka lines emitted in an mCV 
with A'/ wd = 1.0 Mq, rh = 1 gcm~ 2 s _1 , e s = and i = 45°. 
The dotted curves show the corresponding profiles for the 
case where photons are injected at a single height (where the 
emissivity peaks) for the same mCV parameters. These pro- 
files are the same as those shown in Fig. [3^, for the i = 45° 
case (middle panel). 

The line profiles calculated using a realistic emissivity 
profile (solid curves in Fig. [7} show some additional smear- 
ing, particularly redward of the line centre, compared to 
the single injection site case (dotted curves in Fig. [7}. This 
can be attributed to a small fraction of photons now being 
emitted from regions closer to the shock, where the tem- 
perature is « 34keV and the flow speed is ~ 1230 km.s . 
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Thus, these photons are more affected by Doppler effects 
than other photons emitted from regions further away from 
the shock. The small smearing redward of the line centre is 
therefore due to a combination of Doppler shift and ther- 
mal broadening. The effect is not prominent because only 
a small fraction of photons are emitt ed close to the shock, 
as predicted by the emissivity profile (|Wu et al.l l200D. The 
FWHM for the 6.675 keV and 6.97 keV line is « 8 eV and w 
9 eV respectively. This corresponds to an additional broad- 
ening of approximately 25% for the 6.675 keV line and 11% 
for the 6.97 keV line. 



3.5 Comparison with Observations 

Chandra /H ETGS observatio n s of a number of mCVs were 
reported by iHellier fc Mukail (|2004 ). The highest signal-to- 
noise spectrum was obtained in two observations of the mCV 
GK Per during its 2002 outburst. This system is classified 
as an intermediate polar (IP); the WD accretor is though t 
to have a magnetic field of « 1 MG llHellier fc Mukaill2004) 
and a mass of M„ d 2 0.87 ± 0.24 Mq (jMorales Rueda et all 
120021 '). Its typical X -ray luminosity i n the 0.17—15 keV band - 
pass is L x « 7.14 x 10 33 ergs _1 (|Vrielmann et al.ll2005h . 
which places a lower limit on the specific mass accretion 
rate, rh ^ 7.3 g cm" 2 s -1 (adoptin g i? wd « 6.7 x 10 8 cm, in- 
ferred using the lNauenberel (|l972h mass-radius relation, and 
assuming that accretion proc eeds onto a fraction ~ 10~ 3 of 
the WD surface area; see e.g. iFrank. King fc Rainel [20021 ). 

Here we present a pilot study of GK Per and produce a 
simple comparison between the profiles of the simulated and 
observed spectra. A more detailed analysis of the data is left 
for future study. The summed first-order HETGS spectrum 
from the 2002 Chandra observations of GK Per is shown 
in Figure [HI The He-like and fluorescence Fe Kq lines are 
detected, and there is a possible excess at the expected en- 
ergy for the H-like line at 6.97 keV. Also shown in Fig. [S] are 
the upper and lower mass simulated HETGS spectra from 
a 10 ks Chandra observation of an mCV with parameters 
appropriate for GK Per: an upper mass of Af w d = 1.11 Mq 
with _R w d = 5x 10 s cm, a lower mass of M w( j = 0.63 Mq with 
Rvd = 8.2 x 10 s cm and for both masses m = 10gcm~ 2 s" 1 
and e B = 4 x 10 -5 . We adopted inclination angles aver- 
aged over an orbital period of GK Per (50° i ^ 90°) 
(Hclli er. Harmer fc Bendm orc 200_3). We also added a nor- 
malized bremsstrahlun g continuum with an electron tem- 
perature kT c = llkeV (|Vrielmann et al.ll2005T ). 

The measured equivalent widths of the 6.4, 6.675 and 
6.97 keV lines in the CTiandro/HETGS spectra of GK Per 
are 260, 117 and 80 eV, respectively. The relative strength 
of the lines in the simulated spectra is fixed by the number 
of photons used in each simulation and the assumed contin- 
uum flux level. We used 10 8 photons for all three lines, so 
that their relative strengths are comparable. We intend to 
relax this condition in future work. There is a remarkable 
similarity between the observed and the simulated spectra 
for the 6.4 keV line in particular. Realistically, the actual 
strengths of the 6.675 and 6. 97 keV lines dep end on the ion- 
isation structure of the flow l|Wu et al.ll200ll ). since this de- 
termines the emissivity of the lines. The H-like and He-like 
lines are thus potentially powerful diagnostic tools that can 
be used to deduce the physical properties of the post-shock 
accretion column in mCVs. The 6.4 keV line strength, on the 



other hand, depends on the flux of the illuminating X-rays, 
which is largest near the WD surface. Fluorescent lines may 
also be produced in t he surface atmospheric region around 
the accretion column IjHellier fc M ukai 2004), and this may 
contribute to the observed equivalent width of the 6.4 keV 
line. 

The observed fluorescent line i n GK Per exhibit s a red 
wing extending to 6.33 keV which IHellier fc Mukail (|2004f) 
attribute to Doppler shifts. In general we find that bulk 
Doppler shifts have a negligible effect on Fe Kq lines. Our 
simulated spectra show a weaker red wing arising from 
downscatterings near the base of the accretion column. The 
observed spectrum also exhibits a shoulder extending 170 eV 
redwa rd of the 6.4 keV line centre, which IHellier fc Mukail 
(|2004h suggest may be due to Compton downscattering. Our 
simulations confirm that recoil can indeed affect the fluores- 
cent line when it is emitted at the base of the accretion 
column, although in our spectra (Fig. [SJ there is no clear 
evidence that recoil extends down to 170 eV redward of the 
6.4 keV line. 



4 SUMMARY AND CONCLUSION 

We have investigated line distortion and broadening effects 
due to Compton scattering in the accretion column of mCVs 
using a nonlinear Monte Carlo technique that takes into ac- 
count the nonuniform temperature, velocity and density pro- 
files of the post-shock column. Scattered Fe Kq lines were 
simulated for a range of different physical parameters: white 
dwarf mass, accretion rate, magnetic field strength and in- 
clination angle. The photon source regions in the post-shock 
flow were determined from the ionization structure and the 
effects due to the bulk velocity Doppler shift and thermal 
broadening are also considered. 

We find that line profiles are most affected by Compton 
scattering in the cases of low white dwarf mass, high specific 
mass accretion rate, strong cyclotron cooling and oblique in- 
clination angles. Both a lower white dwarf mass (or equiv- 
alently a larger white dwarf radius) and a higher accretion 
rate result in a higher optical depth, hence more pronounced 
scattering effects. Strong cyclotron cooling associated with 
a high magnetic field strength results in a higher density 
in the post-shock flow and hence, a correspondingly larger 
optical depth. Recoil signatures are evident for oblique view- 
ing angles when cyclotron cooling becomes important. These 
are due to scatterings near the base of the flow, where both 
the dynamical and thermal velocities are small. Sharp up- 
scattering features are also seen in the line profiles for large 
viewing angles. These upscattering features are attributed to 
head-on collisions with hot electrons near the shock. Pho- 
tons emitted close to the shock in the accretion column dis- 
play more thermal Doppler broadening. This is generally 
the case for the Fe xxvi 6.97 keV line. The Fe 6.4 keV and 
Fe xxv 6.675 keV lines are generally emitted at or near the 
base of the accretion column respectively, and show a rel- 
atively small amount of thermal Doppler broadening. Bulk 
velocity induced Doppler shifts are negligible compared with 
scattering. 

We also investigated the effects of dispersing the 6.675 
and 6.97 keV line photons using an emissivity calculated 
self-consistently from the ionization structure of the post- 



Compton Scattering of Fe Ka Lines in mCVs 9 



0.0005 
& 0.0004 



a 
= 

o 

u 

■a 



o 
Z 



0.0003 



2 0.0002 



0.0001 




z_ Low Mass Simulated Spectrum 




Observed Spectrum 
6.1 6.2 6.3 6.4 



6.5 6.6 6.7 
Energy (keV) 



6.8 



6.9 



7.0 



7.1 



Figure 8. Fe Kct emission lines for the mCV GK Per detected by Chandra/HIETGS (bottom) and simulated by our model for a lower and 
upper white dwarf mass M w( j = 0.63 Mq (middle) and M w( j = 1.11 Mq (top) both with a specific mass accretion rate m = 10 g cm -2 s _1 . 
The simulated spectra have been convolved with the Chandra /HETGS response function. Cyclotron cooling is negligible and we have 
averaged the simulated lines over all inclination angles in an orbital period. An offset has been added to the simulated plot to allow a 
comparison between the spectra. See text for further details. 



shock column. We find that, compared to injecting the line 
photons at a single height where the emissivity peaks, us- 
ing an emissivity profile does not significantly change the 
overall profiles of the 6.675 and 6.97 keV lines. The result- 
ing line profiles show a small degree of smearing redward 
of the line center due to additional dynamical and thermal 
Doppler effects associated with a small fraction of line pho- 
tons emitted close to the shock. The effect is most important 
for the 6.675 keV line. We estimate a fractional increase in 
the FWHM of no more than » 25%. 

In general, our simulations predict that recoil in scat- 
tering with cool electrons near the base of the column as 
well as upscattering by hot electrons near the shock can im- 
print signatures on the profile of lines emitted near the base 
of the flow. Such Compton signatures may thus be used to 
determine the primary source region of the fluorescent line. 
We predict that when 6.4 keV photons are emitted in the 
dense, cool plasma at the base of the accretion column, they 
can suffer strong downscattering when propagating through 
the electrons streaming downward in the flow. The scatter- 
ing probability is governed by the effective scattering optical 
depth, which increases with the mass accretion rate of the 
system. 

We convolved simulations with the Chandra response 
function and compared them to Chandra /HETGS observa- 
tions of GK Per. Our results indicate that the red wing seen 
in the 6.4 keV line in GK Per could be attributed to Comp- 
ton recoil near the base of the flow. 

Finally we remark that tighter constraints on the dy- 
namics and flow geometry in magnetized accreting com- 
pact objects can be obtained by considering the polariza- 
tion prope rties of the lines (see ISunvaev fe Titarchuk1ll985l ; 
Matt 2004). There is currently considerable effort to develop 
X-ray polarimeters which c an detect de grees of polarization 
of the order of one percent (|Costa et al.ll200l"h . The spectral 
resolution of these detectors should be adequate to search 



for different polarization degrees in emission iron lines. It is 
possible to include a polarization treatment in our calcula- 
tions, but the computational algorithm in our Monte Carlo 
code will need to be revised, which we leave for a future 
study. 
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APPENDIX A: THE WIDTH OF THE 
COMPTON SHOULDER 

Here, we derive an estimate for the width of the prominent 
wings near the base of the line profiles. Let the accretion 
flow be along the z-axis and consider an observer on the 
x — z plane. Let i be the line-of-sight inclination angle of the 
accretion column, and let /3 be the local velocity of the flow, 
normalized to the speed of light c, which can expressed as 



The normalized vector of a scattered photon of energy E' 
propagating in the direction to the observer is 

k = sini x + cosi z . 

Suppose that the normalized vector of the incident photon, 
with energy E, is 

k = sin 8 cos 4> x + sin 8 sin <j) y + cos 8 z . 

Then the change in energy of the photon, after scattering 
with an electron in the flow, due to bulk motion is given by 



EL 

E 



where 7 = (1-/? 2 )~ 1/2 , fi/3 = k-f3, and /x'/3 = k'-fi . In terms 
of the viewing inclination angle and the photon propagation 
vectors, 

E' l+/3cos6» 

E (1 + /3cosi) + ^ c a (1 — sinisin#cos0 — cosicosfl) 



1 + /3(cos 8 — cos i) — A(l — sin i sin 8 cos (f> — cos i cos 8) 



Then for /3 < 1 and E < m c c 2 

EL 
E 

It follows that 
AE 

— — ~ /3(cos 8 — cos i) — A(l — sin i sin 8 cos cf> — cosi cos 8) 
E 

where AE = E' - E and A = E/m c c 2 . 

The maximum energy downshift is caused by the recoil 
process when the photons are scattered by "cold" electrons 
(with (3 ~ 0). This occurs when <j> — n, 8 — n — i, which 



yields 

AE - 2A 

max down 

This result is practically independent of the WD mass 
and the viewing inclination angle. For the 6.4 keV line, 
A£/E| max down » -2.50 x 10" 2 . 

The maximum energy upshift is caused by a Doppler 
shift when the photons are scattered by the fastest available 
downstream electrons (i.e. /3 is no longer a negligible factor). 
The condition for its occurrence can be derived as follows. 
Set 

d_ I AE\ 
M \~E~) 



= 



P(cos8 — cosi) — A(l =p sinisin^ — cosicosS) . 



which gives two conditions for extrema with respect to the 
azimuthal coordinate: <j> = and cj> — it, corresponding re- 
spectively to 

AE 
~~E 

Differentiating the above expression with respect to 8 and 
setting the resulting expression to zero yields the following 
condition for the extrema: 



cos 8 
where 

i = 



Asini 



f3 + A cos i 

The first case {cj> = 0) leads to a maximum energy upshift, 
which requires 



sin U = 



cos 8 



At the viewing inclination angle i = n/2, 



sm8 



cos 8 = 



A 



Hence, the maximum energy upshift is given by 



AE 
E 



sjp + A 2 — A 



If we assume that the maximum (3 takes the free-fall veloc- 
ity at the WD surface, then j3 is simply the reciprocal of the 
square root of the WD radius in the Schwarzschild unit, i.e. 
(3 = v /2GA/ wd /i?w d c 2 . For Af wd = 1.0 Mq,/3= 2.31 x 10~ 2 , 
and for M„ d = 0.5 Mq, j3 = 1.2 3 x 10~ 2 (assum ing the 
Nauenberg mass-radius relation (Nauenberg] Il972l )). this 
gives AE/E\ m ^ up w 1.36 x 10" 2 and 5.06 x 10" 3 respec- 
tively for the 6.4 keV line. 

The Compton shoulder of an Fe line is due to a single 
scattering event. For the 6.4 keV line, the broadening ex- 
tends over 6.24 keV < E < 6.49 keV for a 1.0 M & WD and 
6.24 keV < E < 6.43 keV for a 0.5 Mq WD (omitting ther- 
mal broadening and flow Doppler broadening) , when viewed 
at i = 7r/2. 
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